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Abstract. A simple linear loop is a simple while loop with linear as- 
signments and linear loop guards. If a simple linear loop has only two 
program variables, we give a complete algorithm for computing the set 
of all the inputs on which the loop does not terminate. For the case of 
more program variables, we show that the non-termination set cannot 
be described by Tarski formulae in general. 
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1 Introduction 

Termination of programs is an important property of programs 
and one of the main research topics in the field of program verifica- 
tion. It is well known that the following so-called "uniform halting 
problem" is undecidable in general. 

Using only a finite amount of time, determine whether a given 
program will always finish running or could execute forever. 

However, there are some well known techniques for deciding ter- 
mination of some special kinds of programs. A popular technique is 
to use ranking functions. A ranking function for a loop maps the 
values of the loop variables to a well-founded domain; further, the 
values of the map decrease on each iteration. A linear ranking func- 
tion is a ranking function that is a linear combination of the loop 
variables and constants. Some methods for the synthesis of rank- 
ing functions and some heuristics concerning how to automatically 
generate linear ranking functions for linear programs have been pro- 
posed, for example, in Colon and Sipma [3J, Dams et al. [4] and 
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Podelski and Rybalchenko [6]. Podelski and Rybalchenko [6] pro- 
vided an efficient and complete synthesis method based on linear 
programming to construct linear ranking functions. Chen et al. [2] 
proposed a method to generate nonlinear ranking functions based on 
semi-algebraic system solving. The existence of ranking function is 
only a sufficient condition on the termination of a program. There 
are programs, which terminate, but do not have ranking functions. 
Another popular technique based on well-orders, presented in Lee 
et al. jS], is size-change principle. The well-founded data can ensure 
that there are no infinitely descents, which guarantees termination 
of programs. 

For linear loops, some other methods based on calculating eigen- 
vectors of matrices have been proposed. Tiwari [7] proved that the 
termination problem of a class of linear programs (simple loops 
with linear loop conditions and updates) over the reals is decidable 
through Jordan form and eigenvector computation. Braverman [lj 
proved that it is also decidable over the integers. Xia et al. [8j con- 
sidered the termination problems of simple loops with linear updates 
and polynomial loop conditions, and proved that the termination 
problem of such loops over the integers is undecidable. In [9], Xia et 
al. provided a novel symbolic decision procedure for termination of 
simple linear loops, which is as efficient as the numerical one given 
in m- 

A counter-example to termination is an infinite program exe- 
cution. In program verification, the search for counter-examples to 
termination is as important as the search for proofs of termination. 
In fact, these are the two folds of termination analysis of programs. 
Gupta et al. [10] proposed a method for searching counter-examples 
to termination, which first enumerates lasso-shaped candidate paths 
for counter-examples and proves the feasibility of a given lasso by 
solving the existence of a recurrent set as a template-based constraint 
satisfaction problem. Gulwani et al. [11] proposed a constraint-based 
approach to a wide class of program analyses and weakest precon- 
dition and strongest postcondition inference. The approach can be 
applied to generating most-general counter-examples to termination. 

In this paper, we consider the set of all inputs on which a given 
program does not terminate. The set is called NT throughout the 
paper. For simple linear loops, we are interested in whether the NT 



is decidable and how to compute it if it is decidable. Similar prob- 
lems was also considered in |12] . Our contributions in this paper 
are as follows. First, for homogeneous linear loops (see Section [2] for 
the definition) with only two program variables, we give a complete 
algorithm for computing the NT. For the case of more program vari- 
ables, we show that the NT cannot be described by Tarski formulae 
in general. 

The rest of this paper is organized as follows. Section [5] intro- 
duces some notations and basic results on simple linear loops. Sec- 
tion [3] presents an algorithm for computing the NT of homogeneous 
linear loops with only two program variables. The correctness of the 
algorithm is proved by a series of lemmas. For linear loops with more 
than two program variables, it is proved in Section H] that the NT 
is not a semi-algebraic set in general, i.e., it cannot be described by 
Tarski formulae in general. The paper is concluded in Section |5j 

2 Preliminaries 

In this paper, the domain of inputs of programs is R, the field 
of real numbers. A simple linear loop in general form over R can be 
formulated as 

PI : while (Bx > b) {x := Ax + c} 

where 6, c are real vectors, A nxn , B mxn are real matrices. Bx > b 
is a conjunction of m linear inequalities in x and x := Ax + c is a 
linear assignment on the program variables x. 

Definition 1. [7J The non-termination set of a program is the set 
of all inputs on which the program does not terminate. It is denoted 
by NT in this paper. 

In particular, 

NT(P1) = {x G R n |Pl does not terminate on x} . 
We list some related results in [7j. 
Proposition 1. [7] For a simple linear loop PI, the following is true. 
— The termination of PI is decidable. 



— If A has no positive eigenvalues, the NT is empty. 

— The NT is convex. 

In this paper, only the following homogeneous case is considered. 

P2 : while (Bx > 0) {x := Ax} . 

Let Bi, ... , B m be the rows of B. Consider the following loops 

U : while (BiX > 0) {x := Ax} . 

Obviously, NT(P2)=P|™ 1 NT(Lj). Therefore, without loss of gener- 
ality, we assume throughout this paper that m — 1, i.e., there is only 
one inequality as the loop guard. The following is a simple example 
of such loops. 

while (4„ + , 2 >0) {(:;):= (- 4 ^)(:;)} ■ 
That is B = (4,1), A = 

3 Two-variable case 

To make things clear, we restate the problem for this two- variable 
case as follows. 

For a given homogeneous linear loop P2 with exactly two program 
variables and only one inequality as the loop guard, compute NT(P2). 

For simplicity, we denote the program variables by X\, X2 and use 
NT instead of NT(P2) in this section. If at is a non-zero point in the 
plane, we denote by c£ a ray starting from the origin of plane and 
going through the point a. 

Proposition 2. NT must be one of the following: 

(1) an empty set; 

(2) a ray starting from the origin; 

(3) a sector between two rays starting from the origin. 

Proof. We view an input (xi,x 2 ) as a point in the real plane with 
origin O. If there exists a point M(x\,X2) € NT, any point P on 




the ray OM can be written as P = kM = (kxi, kx2) for a positive 
number k. So BA n (kxi, kx<i) T = k n BA n (x\,X2) T > for any n G 
N. That means P G NT. Therefore, it is clear from the item 3 of 
Proposition [1] that the conclusion is true. 

By the above proposition, the key point for computing the NT is to 
compute the ray(s) which is (are) the boundary of NT. We give the 
following algorithm to compute the ray(s) (and thus the NT) for P2 
if the NT is not empty. The algorithm, as can be expected, is mainly 
based on the computation of eigenvalues and eigenvectors of A. The 
correctness of our algorithm will be proved by a series of lemmas 
following the algorithm. 



Algorithm 1. NonTermination 



Input: Matrices ^2 X 2 and -Bi X 2- 
Output: The NT of P2 with A and B. 

1 if A = or B = then 

2 |^ return 0; 

3 Compute the eigenvalues of A and denote them by Ai, A2; 

4 if Ai -f A A 2 i> then 

5 |^ return 0; // Proposition [T] 

e Take a £l 2 \ {0} such that Ba = and BAa > 0; 
7 if BAoto = then 

choose £ such that B£ > 
if B(A£) > then 

| return {x\x G M. 2 ,Bx > 0} // Lemma Q] 
else 

return // Lemma [5] 



9 
10 
11 
12 



13 if Ai = V A 2 = then 

14 [_ return {x\x G R 2 , Bx > 0, BAx > 0}; // Lemma |6] 

15 Suppose Ai > A 2 

16 if Ai > A2 > then 



17 
18 



choose an eigenvector (5^ related to A 2 such that Bfiz > 0; 

return {x\x = fcicco + fc2/3 2 , fci > 0, > 0}; / / Lemmas [7] and [8] 



19 if Ai > A A 2 < then 
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22 
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24 



if Ai > I As I then 

I let a-i — j4 _1 q:o and return {x\x = fcicco + kia—i, ki > 0, > 0}; 

if Ai < I Aa I then 

choose an eigenvector /3 related to Ai such that B(3 > and 
return {x\x — k/3, k > 0} // Lemma [TUI 



Ax, 

NT 
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Fig. 1. Lemma[T] 



Lemma 1. Suppose NT is not empty and <9NT is the boundary of 
NT. If x G <9NT and Bx ^ 0, t/ien Ax G <9NT. 

Proof. Obviously, S is a linear map from R 2 to R . Because -By > 
for all y G NT, we have Bx > 0. And thus Bx > by the assumption 
that Bx ^ 0. Hence, there exists an open ball Oi(x,ri) such that 
By > for all y G Oi(x,ri). 

Let F be the linear map from R 2 to R 2 that F(y) = Ay for 
any y G R 2 and hence F is continuous. So for any neighborhood 
o(Ax,r) of Aa?, there exists a positive real number r2 such that 
02(x,r 2 ) C oi(a3,ri) and F(o 2 (x, r 2 )) Q o(Ax 1 r). Because a? G <9NT, 
there exist y, z G o 2 {x,r2) such that y G NT and z ^ NT. Then 
A{y), A(z) G o(Ax,r), A(y) G NT and A(z) (£ NT. It is followed 
that there are both terminating and non-terminating inputs in any 
neighborhood of Ax. Therefore, Ax G <9NT. 

Lemma 2. Suppose NT is neither empty nor a ray and <9NT n 
{x\Bx = 0} = {(0,0)}. If By = and BAy > 0, then Ay G NT. 

Proof. By Proposition [2j <9NT consists of two rays. Let h,l 2 be the 
two rays. Since neither l\ nor l 2 is on Bx = 0, l\ and /2 are not 
collinear. So we can choose two points z G l\ and v G l 2 such that 
fiz > 0, Bv > and y = t 1 z + t 2 v for some t x G E, t 2 G E. By 
Lemma (TJ Az and must be on the boundary of NT, i.e., li or l 2 . 
Thus, we have at most four possible cases as follows. 
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Fig. 2. Lemma [2] 



(1) Az = kiZ,Av = k 2 v, (i.e., Az G li,Av G / 2 ) 

(2) = k x z,Av = fc 2 z, (i.e., G G k) 

(3) Az = fcji^AtJ = k 2 v, (i.e., Az G l%,Av G l 2 ) 

(4) Az = kiV,Av = k 2 z, (i.e., Az G Z 2 ,At> G Zi) 

where k\ > 0,k 2 > 0. 

Case (1). Because By = t\Bz + t 2 Bv = and 

BAy = BA(Uz + t 2 «) = t x k x Bz + t 2 /c 2 Pv > 0, 

we have t\t 2 < 0. Without loss of generality, assume that t\ > and 
* 2 < 0. We denote t x Bz by P. Note that P > and t 2 Pv = —P. 
Since PA?/ = (&i — A; 2 )P > 0, we have fci > k 2 > and 

PA n (Ay) = k™ + XBz + k^ + H 2 Bv = k^ +1 P - k r 2 l+1 P > 

for any n G N. By the definition of NT, Ay G NT. 

Case (2). Because BAy = (tiki + t 2 k 2 )Bz > 0, we have 

BA n (Ay) = k^{tiki + t 2 /c 2 )Pz > 

for any n G N. By the definition of NT, we have Ay G NT. 
Case (3). Similarly as Case (2), we can prove Ay G NT. 
Case (4). We shall show that this case cannot happen. Let 

S = {x\x = tiy + r 2 Ay, r\ > 0, r 2 > 0} 



be the sector between the two rays "?/ and Ay. For any w G S, we 
have Sit; = rxfij/ + r 2 BAy = r 2 BAy > 0. 
Because 

A 2 y = A{t x kiv + t 2 k 2 z) = ti~ki~k 2 z + t 2 k 1 k 2 v = hk 2 y, 

we have Aw = r\Ay + r 2 A 2 y = r\Ay + r 2 k\k 2 y G 5*. Therefore, 
G NT and S C NT. As "?/ is a boundary of 5* and i?y = 0, 

is contained in <9NT, which contradicts with the assumption of the 

lemma. So (4) cannot happen. 
In summary, Ay G NT. 




Fig. 3. Lemma 



Lemma 3. If dNT is composed of two rays l\ and l 2 , then either l± 
or l 2 is on Bx = 0. 

Proof. Assume neither l\ nor l 2 is on Bx = 0. Choose a point y such 
that y ^ , By = and BAy > 0. 

Suppose BAy = 0. As NT is not empty, there exists z G NT. 
Hence Ay can be rewritten as Ay = hiz + h 2 y for some hi 6 K, /i 2 6 
R. As a result of BAy = h x Bz + h 2 By = h x Bz = 0, hi = 0. Note 
that 

A n y = h n 2 y,BA n y = h n 2 By = Q . (1) 

According to Eq.flT]) and z G NT, we have BA n (k 1 z + A; 2 y) = 
k x BA n z + k 2 BA n y = k 1 BA n z > for any fci > 0,n 6 E Hence 



{x\x = faz + k 2 y, h > 0} C NT. Therefore, {x\Bx = 0} = <9NT, 
which contradicts with the assumption. 

If B Ay > 0, Ay G NT follows from Lemma H Let S* = {a^iy + 
fc 2 At/, fci > 0, £; 2 > 0}. And we have BA n 2 = k 1 BA n y + k 2 BA ll+1 y > 
for any n G N, z G S. Thus z G NT and S 1 C NT. By the method of 
choosing y, C <9NT. That means is Zi or l 2 , which contradicts 
with the assumption. 

Lemma 4. Suppose A has positive eigenvalues and has an eigen- 
vector a. satisfying Bet = 0. If £ is a vector such that B£ > and 
BA£ > 0, then NT = {x\Bx > 0}. 

Proof. For any y G {x\Bx > 0}, it can be written as y = k\£ + k 2 oc 
for some k x G R, k 2 G R. As By = k x B^ + k 2 Ba = k x B^ > 0, we 
have ki > 0. Thus B Ay = k 1 BA^ + k 2 BAoL = k x BAi > and Ay G 
{x\Bx > 0}. By the definition of NT, we have {x\Bx > 0} C NT 
and hence NT = {x\Bx > 0}. 

Lemma 5. Suppose A has positive eigenvalues and has an eigen- 
vector a satisfying Bol = 0. If there is a vector £ such that B£ > 
and BA£ < 0, then NT = 0. 

Proof. For any y G {x\Bx > 0}, it can be written as y = k x a + k 2 £ 
for some k\ G R, k 2 G R. Since By = k 2 B^ > 0, we have k 2 > 0. 
And because BAy = k 2 BAi < 0, NT = 0. 

Lemma 6. Suppose A has a positive eigenvalue and a zero eigen- 
value. If 7 is an eigenvector related to the positive eigenvalue such 
that Bj > 0, then NT = {x\Bx > 0, BAx > 0}. 

Proof. Let (3 be an eigenvector with respect to eigenvalue and A be 
the positive eigenvalue. Let S be the set {x\Bx > 0, BAx > 0}. For 
any y G S, it can be written as ki(3+k 2 ~f for some k\ G R, k 2 G M. We 
have BAy = k 2 XBj > 0, thus k 2 > 0. Note that BA n y = k 2 \ n j > 
for any n G N, hence S C NT. Because {cc|i?a; < V BAx < 
0} n NT = 0, NT = {x\Bx > 0, BAx > 0}. 

Lemma 7. Suppose A has two positive eigenvalues Ai > X 2 > and 
two eigenvectors /3i and ^2 related to Ai and A 2 , respectively, such 
that B(3i > 0,5/32 > 0. If ex. is a vector such that Bex = and 
BAa. > 0, then NT = {x\x = hot + k 2 f3 2 , h >0,k 2 > 0}. 



Proof. It is easy to know (3i,(3 2 G NT, thus NT is neither empty 
nor a ray. By Lemma [3] there is a "?/ C <9NT and y satisfies By = 0. 
Since for any z G <9NT, we have BAz > 0. So BAy > and hence 
c£ = In other word, c£ is one ray of <9NT. Let the other ray 
of <9NT be I. As — BAck < 0, —a is not Z. By Lemma [H we have 
Al G <9NT. So Z is one of ^i,(3 2 and A~ x a. By directly checking, we 
know j3 2 is Z and so NT = {x\x = fe x a + fc 2 /32> &i > 0, Zc 2 > 0}. 

Lemma 8. Assume that A has one positive eigenvalue A with mul- 
tiplicity 2 and only one eigenvector (3 satisfying B(3 > 0. If ct is 
a vector such that Bex = and BAa > 0, then NT = {x\x = 
hio. + h 2 (3, ki >0,k 2 > 0}. 

Proof. By the theory of Jordan normal form in linear algebra, there 
exists a vector (3\ such that Af3\ — (3 + A/3i and (3 and /3i are 
linearly independent. 

Let cki = Aol. We claim that 

\/n G n.{BA n otx > A 3h 2 > 0.(A n ai = h/3 + h 2 (3 1 )). (2) 

To prove this claim we use induction on the value of n. 

Suppose a = h\f3 + h 2 (3 1 . If n = 0, then a± = Aa = (Z^A + 
h 2 )(3 + h 2 \\3\. Because Ba x = XBa + h 2 B(3 = h 2 B(3 > 0, we have 
h 2 > 0. 

Now assume that the claim is true for n — 1. Let A n ~ 1 a\ = h\(3 + 
h 2 (3i where h 2 > 0. Because A n ax = A(A n ~ 1 a 1 ) = (Xhi + h 2 )f3 + 
AZi 2 /3i, we have \h 2 > and BA n a 1 = \BA n ~ 1 a 1 + h 2 B(3 > 0. So 
the claim is true for any n G N and we have a^ G NT. 

Obviously, /3 G NT and /3 and ai are linearly independent, so 
NT is not a ray. By Lemma El c£ C <9NT. 

Let the other ray of <9NT be Z. As —BAa < 0, — o : is not Z. By 

Lemma [U Al = I or Al = c£. So Z must be ~p or A -1 ai. By directly 
checking, we know Z is p and thus NT = {x\x = kict + k 2 (3, kx > 
0,fc 2 > 0}. 

Lemma 9. Suppose A has a positive eigenvalue Ai and a negative 
eigenvalue X 2 with Ai > |A 2 | and two eigenvectors (3\ and (3 2 related 
to Ai and \ 2 , respectively, such that B(3\ > 0,Bf3 2 > 0. Suppose 
a is a vector such that Bex = and BAa > 0. Let a_x = A~ 1 a, 
e*i = Aa. Then NT = {k ± a + fc 2 a-i, fci > 0, k 2 > 0}. 



Proof. Let ck_i = h±(3i + hifii- So a = Aa_i = h\\\fi\ + h 2 \ 2 (3 2 
and Qi = Aol = hi\\(3i + /i 2 A|/3 2 - Because Bq = and Ba 1 > 0, 
/ii, h 2 and Aa_! are all positive. 

Note that aix = (— AiA 2 )a_i + (Ai + A 2 )a: where — AiA 2 > and 
Ai + A 2 > 0. Let S = {x\x = kict + fc 2 a_i, k\ > 0, k 2 > 0}. Since 
By = k 2 Bct_ 1 > and Ay = (k 2 + fci(Ai + A 2 ))a — A;iAiA 2 a_i G S 
for any y G S, we have NT D 5. 

Let y = kia+k 2 OL- 1 . Because By = k 2 Ba_ 1 < for any k 2 < 
and BAy = k 1 Ba 1 < for any k x < 0, we have NT = S. 

Lemma 10. Suppose A has a positive eigenvalue \± and a negative 
eigenvalue A 2 such that Ai < |A 2 |. If there are two eigenvectors (3\ 
and 02 related to Ai and X 2 , respectively, such that B(3± > and 
Bf3 2 > 0, then NT = {x\x = k/3 1: k > 0}. 

Proof. Consider any (3 = k\ft\ + A; 2 /3 2 G M 2 . 

If k 2 ^ 0, because A n (k 1 (3 1 + k 2 /3 2 ) = fciAJ^i + fc 2 A^ 2 and 

BA n (k 1 f3 1 + k 2 p 2 )BA n + x {k x f3 x + A; 2 /3 2 ) < 

when n is large enough, + A; 2 /32 ^ NT. 

If k 2 = 0, obviously, NT D {x\x = kfa, k > 0} and Bkfo <£ NT 
for any k < 0. 

So NT = {x|£c = k(3 1 ,k> 0}. 

Now, the correctness of our algorithm NonTermination can be 
easily obtained as follows. 

Theorem 1. The algorithm NonTermination is correct. 

Proof. First, the termination of NonTermination is obvious because 
there are no loops and no iterations in it. Second, it is also clear that 
the algorithm discusses all the cases of eigenvalues of A, respectively. 
According to Lemmas 4-10 (each of them corresponds to a certain 
case in the algorithm as commented in the algorithm), the output of 
the algorithm in each case is correct. 

Example 1. Compute the NT of the following loop. 
w M e(4, 1+ , 2 >0) {(£) = 

Herein, B = (4, 1),A = 




The computation of NonTermination on the loop is: 
Line 1. B ^ and A ^ 0. 
Line 4. A has a positive eigenvalue — 1 + \/l7. 
Line 6. Let a = (-1, 4) T , (*i = Aa = (18, -4) T 
Line 7. 5c*! = 68 ^ 0. 

Line 13. The two eigenvalues of A are — 1 + y/V7,—l — ^/Tf, 
respectively. Neither of them is 0. 

Line 19. A has two eigenvalues, of which one is positive and the 
other negative. 

Line 20. The absolute value of the negative eigenvalue is greater 
than the positive eigenvalue. 

Line 22. The eigenvector with respect to the positive eigenvalue 
is (3 = (1, v5J+!) T anc i B p > . Return {x\x = k(3,k> 0}. 

4 More variables 

Theorem 2. In general, NT is not a semi- algebraic set. 

Remark 1. All Tarski formulae are in the form of conjunctions or/and 
disjunctions of polynomial equalities and/or inequalities, so, in other 
words, semi-algebraic sets are exactly the sets defined by Tarski for- 
mulae. By Theorem El we can conclude that the non-termination 
sets of linear loops with more than two variables cannot be defined 
by Tarski formulae in general. 

Remark 2. It should be noticed that all polynomial invariants are 
semi-algebraic sets. 

In order to prove the above theorem, we give an example to 
demonstrate its NT is not a semi-algebraic set. 

Proposition 3. Let a linear loop with three program variables be as 
follows. 



P3 : while (x\ + 2x 2 + x 3 > 0) 




Then NT(P3) is not a semi-algebraic set. 



The conclusion can be proved by using the following lemmas. For 
simplicity, NT(P3) is denoted by NT in this section. 

Lemma 11. Denote by r the following set 

{9(xl + x\) - x\ < 0, x 3 > 0}, 

then t C NT. 

Proof. For any (xi,x 2 ,x 3 ) G r, we have x 3 > 3|xi|,£3 > 3|x2| and 
thus X\ + 2x 2 + x 3 > 0. Because A(xi, x 2 , x 3 ) T = (2xi, 3x 2 , 5x 3 ) T and 
9(4x? + 9x 2 2 ) ~ 25x l < °> A(x u x 2 ,x 3 ) T G r. Therefore r C NT. 

Lemma 12. <9NT C NT. 

Proof. Because the loop guard is of the form B(x\, x 2 , x 3 ) T > 0, NT 
is a closed set. So the conclusion is correct. Furthermore, for any 
(xi, x 2 , x 3 ) G <9NT, x\ + 2x 2 + x 3 > 0. 

Lemma 13. If (x\,x 2 ,x 3 ) G NT and A(x\, x 2 , x 3 ) T G <9NT, then 
(xi,x 2 ,x 3 ) G <9NT. 

Proof. Let a; = (xi,x 2 ,x 3 ). If the conclusion is not true, there exists 
a ball o(x,r) C NT. Because v4a? T G <9NT, there exists a;' such that 
\Ax — x'\ < r and x' is not in NT. 

Since \A~ 1 x' — x\ < \x' — Ax\ < r, A~ 1 x' G o(x, r). So A~ l x' G 
NT and thus x' G NT, which is a contradiction. 

Lemma 14. {(i -i C «9NT. 

Proof. Let p n = —3^, ^r),w > 0. We use induction on the value 
of n. 

When n = 0, because -£>p — -B(l> — 1> 1) T = and 

SA fc p = 2 fc - 2 x 3 fc + 5 fc > for any jfc G N + , 

we have po £ <9NT. 

Now assume that the conclusion holds for n— 1. So, v4p ra = p n -i G 
<9NT C NT. By Lemma O p n G <9NT. 

Lemma 15. For any non-zero polynomial f(xi,x 2l x 3 ) G M[xi, x 2 , x 3 ], 
there exists an N such that /(^r, — ^) 7^ /or a// n > N. 



Proof. Assume that the conclusion does not hold. Then there exists 
a subsequence {((|) nfc , — (|) nfc , (l) Uk )}T=i sucn that / vanishes on 
each point of it. 

Let / = bxx^x^xf + ... + b s x1 a x p 2 8 xl s where b { G MA ^ 0, a, t G 
N, 0i G N, 7i G N, and (a*, A, 7i) 7^ (ay, ft, 7i) for « ^ J- 

Obviously s > 1 because / ^ 0. Let = (|) a< (§) ft (J) 7i - 

It is an obvious fact that 2^3^ 5^ ^ 2 a *3 ft 5 7 * for i ^ j. Hence 
ti,t 2 , are pairwise distinct. Without loss of generality, let ti > 
t 2 > ... > t s . 

For every j > 1, we have lim (j-) nk = 0. Thus 

fc— >oo 1 
f((i\ n k —(l)n k (l\n k \ 

This contradicts with f{{\) nk , -(|) nfc , (^) nfc ) = 0. Therefore the 
conclusion follows. 

Using the above lemmas, we can now prove Theorem [2j 

Proof. Denote by S the sequence {(^) n , — (|) n , (|) n )}- By Lemma 
MS C <9NT. 

Assume NT is a semi- algebraic set. Then there exist finite many 
polynomials /jj G M[xi, x 2 , x 3 ] and < i(J - G {<, =} for z = 1, s and 
j = 1, rj such that 

S Ti 

NT = Uf|{(x 1 ,x 2 ,x 3 )Gf 3 |/ !J ^0}. (3) 
i=ij=i 

Because S C <9NT C = 0}jj, for any x E S, there exists a 
polynomial /y such that fi,j(x) = 0. By pigeonhole principle there 
exists an /y and a subsequence Si of S such that /jj vanishes on 
Si, which contradicts with Lemma [151 

5 Conclusion 

In this paper, we consider whether the NT of a simple linear 
loop is decidable and how to compute it if it is decidable. For ho- 
mogeneous linear loops with only two program variables, we give a 
complete algorithm for computing the NT. For the case of more pro- 
gram variables, we show that the NT cannot be described by Tarski 
formulae in general. 
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